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<n: We study hidden parasupersymmetry structures in purely bosonic quantum mechanics on 

__■ compact equilateral graphs. We consider a single free spinless particle on the graphs and show 

that the Rubakov-Spiridonov-Khare-Tomiya parasupersymmetries are hidden behind degener- 
"^T* 1 ate spectra. 
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1 Introduction 



Parasupersymmetry is a generalization of supersymmetry to parastatistics: It is a symmetry between 
bosons and parafermions. Historically, parasupersymmetry was first introduced in non-relativistic 
quantum mechanics by Rubakov and Spiridonov [1], and nowadays there are two different for- 
mulations of parasupersymmetric quantum mechanics. The one is parasupersymmetric quantum 
mechanics of order 2 proposed by Rubakov and Spiridonov [1], which is extended to arbitrary order 
by Khare [2] and Tomiya [3], and the other is that proposed by Beckers and Debergh [4], which is 
extended to arbitrary order by Chenaghlou and Fakhri [5]. Common ingredients of these two parasu- 
persymmetric quantum mechanics are Hamiltonian operator H, parafermion number operator Npf, 
nilpotent parafermionic charge (parasupercharge) Q + and its adjoint Q that satisfy (Q ± ) p+1 = 0, 
where p G N is the order of parafermion. Parasupersymmetric quantum mechanics is characterized 
by the parasupersymmetry algebra that consists of the commutation relations 



and some multilinear relations among H, Q + and Q , which depend on the formulations but are 
reduced to the ordinary jY = 2 supersymmetry algebra when p = 1. The most important feature 
of parasupersymmetric quantum mechanics of order p is that the energy spectrum exhibits (p + 1)- 
fold degeneracy. This can be easily seen from the commutation relation (1.1b) and the nilpotency 
(Q ± ) p+1 = 0: Since Q + and Q raise and lower the parafermion number by one, by defining a 
boson state as an annihilated state by Q , one immediately sees that successive multiplications of 
Q + to the boson state yield p distinct parafermion states with the same energy eigenvalue. 

The purpose of this paper is to explore parasupersymmetry structures hidden behind degenerate 
spectra in single-particle quantum mechanics on graphs (quantum graphs). In order to realize para- 
supersymmetric quantum mechanics of order p, we need to have some symmetry transformation 
which plays the role of parafermion number operator, or grading operator. Such symmetry transfor- 
mation must commute with the Hamiltonian and have p + 1 distinct eigenvalues, because Npf splits 
the total Hilbert space into p + 1 distinct subspaces. A typical example of such symmetry transfor- 
mations on graphs is Z p+ i cyclic rotation of p+ 1 edges, Z : y(xj) \— > y{ x j+i) (mod p+l), whose 
eigenvalues are {l,q, • • • ,q p } with q being the (p+ l)fh root of unity. (Here y(xj) is a wavefunction 
on the jth edge.) If the Hamiltonian H is independent of the edges, H and Z obviously commute. 
Thus the cyclic rotation Z is a good candidate for the grading operator q NpF , up to the question 
whether the boundary conditions are invariant under Z. Notice that the cyclic rotation Z = ^ A ' PF is a 
natural generalization of ordinary parity operator V : x \- > —x, which plays the role of fermion parity 
(— l) F in one-dimensional quantum mechanics on M with hidden ,JV = 2 supersymmetry [6]. Note 
also that the cyclic rotation of edges is well-defined only for equilateral graphs. 

The rest of this paper is organized as follows. In Section 2 we set up our model. In this paper 
we consider compact equilateral graphs that consist of two vertices connected by N edges of equal 
finite length. For the sake of simplicity in this paper we focus on the free particle Hamiltonian (or 
minus the Laplace operator), which is obviously edge independent. We then impose a generic Zjv 
cyclic symmetry in this system and obtain a Z^-graded Hilbert space. In Section 3. 1 we show that 
the Rubakov-Spiridonov-Khare-Tomiya (RSKT) parasupersymmetry algebra of order p = N — 1 
is hidden behind the spectral problem with A^-fold degeneracy. We classify boundary conditions 
invariant under parasupersymmetry transformations in Section 3.2 and then explicitly derive the 
parasuperspectrum in Section 4. It turns out that there are only four types of parasupersymmetry 
invariant boundary conditions, of which two leads to spontaneous breaking of parasupersymmetry. 
Section 5 is devoted to conclusions and discussions. 



[H,Q ± ] = [H,N VP }=0, 

[^vp F ,e ± ] = ±e ± , 



(l.la) 
(1.1b) 
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(a) Unfolding picture. (b) Folding picture. 

Figure 1: (a) Carambola graph Cn and (b) its folding. Arrows indicate the direction of coordinates. 

2 Z^-graded Hilbert space 

In this section we study non-relativistic quantum mechanics for a free particle on compact graphs 
with Z# symmetry and its grading structure. As mentioned in the previous section, in this paper 
we will focus on equilateral graphs that consist of two vertices connected by N edges of the same 
length i. Following Ref. [7] we will call this graph a carambola graph and denote it by C^} An 
important point to note is that, by folding the edges, quantum mechanical system on such graphs 
has an alternative equivalent description as a vector quantum mechanics on a single interval (0,£); 
see Figure 1 . This is referred to as the folding trick in the boundary conformal field theory literature 
[8,9]. As shown in Figure 1(b), wavefunction of the system is given by an Af-component complex 
column vector in the folding picture, \ff(x) = (Yi (x), ■ ■ • , Yn(x)) t , where T stands for transposition. 
Namely, the Hilbert space is given by a tensor product 

H=L 2 (0,£)®C N . (2.1) 

A free particle on the graph is described by the time-independent Schrodinger equation H\jr(x) = 
Ey(x), where H is the Hamiltonian operator given by the N x N diagonal matrix 

d 2 

H = --^®I, (2.2) 

with / being the N x N identity matrix. (Throughout this paper we will work in the units h = 
2m = 1.) Given an observable, we have to specify the domain on which it becomes the self-adjoint 
operator. For the case of the Hamiltonian, this can be done by specifying boundary conditions for 
iff consistent with the probability current conservations at each boundary (or Kirchhoff 's law of 
probability current at the vertices). It is well-known that probability current conservations leads to 
the following U (N) x U (N) family of boundary conditions (see e.g. Ref. [10]): 

(7 + ^)^(0) + /M (/-£/)y>i;0) =0, UeU(N), (2.3a) 
(I + U)f(£)-iM (I -&)$(£) =0, UeU(N), (2.3b) 

where prime (') indicates the derivative with respect to x. Mq is an arbitrary scale parameter with 
length dimension —1, which must be introduced to adjust the scaling dimensions of the first and 
second terms in the boundary conditions. We note that without any loss of generality Mo can be 
taken to be positive, Mq > 0. The sign difference in the second terms of (2.3a) and (2.3b) is just for 
later convenience. 

'Carambola (also known as star fruit) is a golden-yellow tropical fruit with a star-shaped cross section. 
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Let us next consider a symmetry transformation of edges that plays the role of the grading 
operator ( e 2m /N)NpF_ ^ s ment i one( j m t h e previous section, a typical example of such symmetry 
transformations is a cyclic rotation of edges, y/j ^ V;'+i (mod N). However, it is not necessarily 
to be the cyclic rotation. In this paper we consider a generic Z# symmetry transformation given by 
the following unitary transformation: 

Z : \jr(x) i-> (Z\jr)(x) = Z\jr(x), (2.4) 

where Z£U (N) is an N x N traceless unitary matrix satisfying Z N = I andZ" ^I(n = 1, •• • ,N— 1). 
Z has the eigenvalues {l,g, • • • ,<? A '~ 1 } with q = exp(27li/N) being the Mh root of unity such that it 
has the following spectral decomposition: 

N-l 

Z = £ <fPn, (2-5) 

n=0 

where P n is a hermitian projection operator onto the nth eigenspace. It should be noted that, in terms 
of the unitary matrix Z, such projection operator can be constructed as follows: 

p _ I + q-"Z + q- 2n Z 2 + ■■■+ q-( N - l >Z N - y 

"n — s (2-D) 

It is not difficult to show that thus constructed projection operator indeed satisfies the orthonormality 

relation P n P m = 8 nm P m , completeness relation Pq H h-fV-i = L and hermiticity P^ = P n , where 

the first two relations follow from the identity 1 + q H h q N ~ l = 0. 

Next we impose the system to be invariant (or symmetric) under the transformation Z. To this 
end, it is worthwhile to remember first the meaning of symmetry in quantum graphs, by follow- 
ing the discussion given by Fiilop et al. [11] for the case of U (2) family of point interactions on 
S 1 . The system is said to be symmetric under Z if the transformed state Z\jf satisfies i) the same 
Schrodinger equation and ii) the same boundary conditions as the original ones. The first point 
is trivially satisfied: If iff satisfies the Schrodinger equation Hy = Ey, then the transformed state 
Z\fir also satisfies the same Schrodinger equation H{Zy) = E{Z\jf) because H = —d 2 /dx 2 ®I and 
Z = 1 <g>Z trivially commutes, [H,Z] = 0. The second point is, however, nontrivial: If y satisfies 
(2.3a) and (2.3b), then the boundary conditions for Zlff become 

(I + ZUZ- l )(Zy)'(0) + iM Q (I-ZUZ- l )(Z\(r)(0) = 0, (2.7a) 
(I + ZUZ- l )(Zy)'(l)-iMQ(I-ZUZ- l )(Zy)(e) =0, (2.7b) 

which follow from multiplying Z to Eqs. (2.3a) and (2.3b) from the left. Comparing these to the 
original boundary conditions (2.3a) and (2.3b), we see that Z induces the maps 

f/4z(/Z _1 and U^ZUZ^ 1 . (2.8) 

If U = ZUZ 1 and U = ZUZ~ l (i.e., [Z,U] = [Z,U] = 0), the boundary conditions are invariant 
under Z. To summarize, if Z is a symmetry of the system, the unitary matrices Z, U and U must be 
all simultaneously diagonalizable; in other words, they have to share the same projection operators 
P n . Such unitary matrices enjoy the following parametric representations: 

Af-l 



U = £ exp(j0„)P„, (2.9a) 

Af-1 

^=£exp(ift,)P„, (2.9b) 



n=0 
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where B n and Q n (n = 0, 1, • ■ ■ ,N — 1) are 2N independent parameters. Eqs. (2.9a) and (2.9b) are the 
general solutions to the conditions [Z,U] = [Z,U] = 0. 

Now, both the Hamiltonian and boundary conditions commute with Z, the Hilbert space com- 
pletely splits into N orthogonal subspaces; that is, % can be written as the following direct sum: 

U = U ®Ui®---®U N -i. (2.10) 

where Ti n = {y„ = P n \jr \ iff GH}. Notice that y n G % n satisfies the eigenvalue equation Zy n = 
q n y n . Let us next investigate the boundary conditions for Yn- To this end, let e n be a simultaneous 
eigenvector of Z, U and U that satisfy the orthonormality and completeness relations, e\fi m = 8 nm 

and e§e\ H heAr-ie#-i = /. The projection operator P n can then be written as P n = e n e] x such that 

\fr n = P n \fr is cast into the following form: 

f n (x) = Y n (x)e n . (2.11) 

where Xj/ n = ej t \fr. In what follows we will call the set of complete orthonormal eigenvectors 
{eo, ■ ■ ■ ,?n-i} the (ordered) basis and Xj/ n (x) the component. By multiplying el from the left, the 
boundary conditions (2.3a) and (2.3b) boil down to the following 2N independent Robin boundary 
conditions for the component y n (n = 0, 1 , • • • ,N —\): 

cos <(0) + Mo sin ¥n (0) = 0, (2.12a) 

cos xi/ n (£) -M sin ¥n (£) = 0, (2.12b) 

which follow from the eigenvalue equations for e\, e]JJ = e lQ "e\ and e^U = e ,e "ej. 

To summarize, by imposing the Z/v-symmetry (2.4), we have obtained the Z^-graded Hilbert 
space (2.10) characterized by the eigenvalue q n = ( e 2jri / A ') n of the unitary matrix Z. In parasu- 
persymmetric quantum mechanics language, this eigenvalue corresponds to the exponential of the 
parafermion number: Z counts the parafermion number modulo N. Each subspace H n carries the 
parafermion number n, and specified by the Robin boundary conditions (2.12a) and (2.12b). 

A few more comments are in order about two other symmetry transformations on the graph C^. 

• The first symmetry transformation is a reflection 1Z around the midpoint x = £/2, 72. : jc i— > 
I — x, which just flips the coordinate direction. The action of 1Z on \jf e % is defined as 

K:f(x)^{Kf)(x) = f(£-x). (2.13) 

Obviously this Z2-transformation 1Z commutes with the free Hamiltonian H such that it pre- 
serves the energy spectrum. However, 1Z does not preserve the domain of H. To see this, 
suppose that \jf fulfills the boundary conditions. Then, the transformed state TZrjf satisfies 

-(I + U)(K\pry(£) + iMo(I-U)(Ky)(l)=0, (2.14a) 
-(I + U)(Kf)\0)-iM (I-U)(Tlf)(0) = 0, (2.14b) 

where we have used (TZ\fr)'(x) = — \f/(£ — x). Thus 1Z induces the maps 

U&U and U&U. (2.15) 

lfU = U, the system is said to be 7£-invariant. If two systems are related by the maps, we say 
such systems are 7£-dual. As mentioned above, 1Z is the spectrum-preserving transformation 
such that IZ-dual systems are isospectral. We will encounter this duality in Section 4.3. 
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Given the Z^r-transformation Z, we can introduce another Z/v-transformation X : y(x) i-> 
(X\}r)(x) that satisfies the ^-commutation relation 

ZX = qXZ. (2.16) 

Indeed, X can be defined as (Xy)(x) = Xy(x) with X being an N x N traceless unitary 

matrix given by X = eoe N _ x +e\e*Q + e2e[ H h?w-i?^_ 2 - Thus constructed unitary matrix 

satisfies X N = I and ^-commutation relation ZX = qXZ, so does X. In this paper we will call 
a pair of unitary matrices (X,Z) the Weyl pair, 2 which in the basis {eo, ■ ■ • takes the 

following standard forms: 



(2.17) 





(° 


• 


• 


1\ 









• 


• \ 




1 


• 


• 










q 


• 





X = 





1 • 


• 





, z = 








q 2 ■ 







v> 


• 


• 1 


«/ 









• 


• q N - 1 ) 



By construction it is obvious that X and X^ raise and lower the parafermion number n 

Xe n = e n+ i and X^e n = e n -\ (modAf)- (2.18) 



As we will see in Section 3.1, the unitary matrix X enables us to construct the parasuper- 
charges : 7-L n — > % n ±\ in the basis independent way. 



3 jV = 2 parasupersymmetry of order p = N—l 

Roughly speaking, parasupersymmetric quantum mechanics of order p is just a set of p ordinary 
jV = 2 supersymmetric quantum mechanical systems glued together in an intertwined way. The 
point is a hierarchy of isospectral Hamiltonians based on the factorization method. In this section 
we first try to construct the system to be A^-fold degenerate by using the factorization method, and 
then reveal the underlying parasupersymmetry structure. 

There are two important observations here. The first is that, except for the case d n ,9 n ^ K 
(mod 2n), the boundary conditions (2.12a) and (2.12b) can be compactly written as (A J y M )(0) = 
= (Ag lff n )(£), where Ag are one-parameter family of first-order differential operators defined by 

A^ = ±— +M(6) with M(6) =M tan( - ) . (3.1) 
dx \2/ 

Notice that A e = - A+ e . Note also that A g is a formal adjoint of Ag . 

The second is that the free Hamiltonian —d 2 /dx 2 can be factorized as A^Aq —M 2 (d) for any 
8, which is the heart of hidden parasupersymmetry structures in our model. To be more specific, let 
H n = —d 2 /dx 2 be the Hamiltonian operator for the component \\f n that satisfies H n \\r n = E\\f n . Then 
the above factorization enables us to construct the following hierarchy of Hamiltonians: 

#o =A-A+ -M 2 (a { ) 

Hi =A+A~ -M 2 (a Y ) = A- A+ -M 2 (a 2 ) 

Hi = =A+ 2 A ai -M 2 (a 2 )=A a3 A+ 3 -M 2 (a 3 ) 



Hn-i = =K N _ l K N - l -M 2 (0W-i). 



"In the literature [12] a pair (X* ,Z) is called the Weyl pair. 
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We emphasize that the parameters {ct n } n =i are independent of the boundary condition parameters 
{B n , 9„}^~q at this stage. Parasupersymmetry relations corresponding to this hierarchy of Hamilto- 
nians are as follows: 

A t Yo (x) = JE + M 2 (a 1 )xi/ 1 (x) , 

< _ V (3.2a) 

A ai Yi( x ) = y£' + M 2 (ai)w ) (x), 

A+ yr^x) = JE + M 2 (cc 2 )y 2 (x), 

< _ i (3.2b) 

K 2 Wi{x) = ^E + M 2 {a 2 )yi(x), 



A a N ,W-2W = jE+M 2 (a N - l )\i/N-i(x), 
< i (3.2c) 

A a N -,WN-i{x) = ^E + M 2 (a N - l )\i/N-2(x). 

Now our task is to tune the parameters {6„, 6 n } to be consistent with Eqs. (3.2a)-(3.2c). For exam- 
ple, if we impose the Robin boundary condition (A+y/b)(0) = 0, which corresponds to the choice 
do = CCi ^ % (mod 2k), the parasupersymmetry relation (3.2a) implies Yi must obey the Dirichlet 
boundary condition yi(0) = 0, which corresponds to the choice 6i = % (mod 2%). It follows from 
the relation (3.2b) that the Dirichlet boundary condition ^i(O) = leads to the Robin boundary 
condition for X(f 2 , (A^ 2 V^)(0) = 0, which corresponds to the choice B 2 = — <x 2 / % (mod 2k). This 
procedure can be easily extended for all \\f n . Before completing this program, however, let us check 
the hidden parasupersymmetry at the algebra level. 



3.1 Parasupersymmetry algebra 



■ 4T4± 



The above constructed hierarchy of Hamiltonians based on the factorization —d 2 /dx 2 
M 2 (a) is well-described by the RSKT parasupersymmetry algebra of order p = N — 1. To see 
this, let us first work in the basis {eo,e\,- • • ,e^-i} in which both H and Z become diagonal, H = 



diag(# ,#i, " ,H N -i) andZ = diag(l,<?,- 
as the following standard forms 

/ • 
A 



JV- 1 N 



In this basis we define parasupercharges Q 









A a 2 



\ 









A + 








/o 


A ai 


• 





\ 








A a 2 ' 













• 


• A a N . 


1 


\o 





• 





/ 



(3.3) 



It is straightforward to show that a set of operators {H,Z,Q ± } satisfies the following ,JV 
parasupersymmetry algebra of order p = N — 1 : 

Z N = Id ^ - 



z 1 = z~ 

= 0, 



(<2 d 

ZQ ± 

[H,Q ± ] = [H,Z}=0, 

N f(Q + ) N - l - m (Q-) N -HQ" 

m=0 



N-l 

= Yl[H + M 2 (a n )], 

n=l 



(3.4a) 
(3.4b) 
(3.4c) 
(3.4d) 

(3.4e) 
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where Id is the identity operator on %. Notice that the multilinear relation (3.4e), which is different 
from the standard one, is an alternative characteristic relation of RSKT parasupersymmetry algebra 
proposed by Huang and Su [13]. Note also that the case N = 2 reduces to the ordinary jY = 2 
supersymmetry algebra, Q + Q + Q Q + = H + M 2 {a{), where the origin of energy is shifted by 
M 2 («i). 

Let us next move on to an arbitrary basis. Although the algebra (3.4a)-(3.4e) is independent 
of the choice of the basis, it is more desirable to construct the parasupercharges Q ± in the basis 
independent way. This can done by using the projection operator P n and the unitary matrix X as 
follows: 

N-2 

e + = £A+ +1 8)XP„, (3.5a) 

n=0 
N-l 

QT= £A" ®X*P n . (3.5b) 

n=\ 

Interpretations of these expressions are obvious: When <2 ± act on w = Wo + Wi H h Wn-i S 'H, 

the projection operator P n picks up the mode Wn = \\f n e n . Then X raises the basis vector e n to e n+ \ 
(see Eq. (2.18)) and A„ the component Wn to Y„+\. Similarly, X^ lowers the basis vector e n to 
e„_ i and A+ the component \\f n to Wn- 1 ■ Thus we obtain the following complex (known as the 
"paracomplex" [14, 15]): 

O^Uo^ni ^± ■■■ ^ H N -i 0, (3.6) 

Q e: gj e w _, 

where = Q ± P„. Recall that the projection operators are written in terms of Z; see Eq. (2.6). 
Hence, in order to construct the parasupercharges, we need only the Weyl pair (X,Z) and the first- 
order differential operators . 

3.2 Parasupersymmetry invariant boundary conditions 

The boundary conditions are said to be parasupersymmetry invariant if the transformed states Q y n 
satisfy the same boundary conditions as those for y n ±i- In this section we classify such parasuper- 
symmetry invariant boundary conditions in our model. 

To this end, we first note that, once we fix the boundary conditions for the bosonic sector 
Mo 3 Wo, the boundary conditions for the parafermionic sectors 1~L n (n = l,--- ,iV— 1) are auto- 
matically determined via parasupersymmetry relations (3.2a)-(3.2c). Note also that there are two 
possibilities for the bosonic sector: the Robin boundary condition (A+ yb)(0) = or Dirichlet 
boundary condition i//b(0) = 0, the former corresponds to the choice Bq = a,\ ^ % (mod 2%) and 
the latter do = 7t (mod 2%). Thus there are two distinct types of boundary conditions at each bound- 
ary, one is the sequence of boundary conditions that starts from the Robin boundary condition for 
Wo and the other the Dirichlet boundary condition. Let us first study the former type. 

Type R. If we impose the Robin boundary condition (A+ Wo) (0) = 0, then the parasupersymmetry 
relation (3.2a) implies that G Hi must satisfies the Dirichlet boundary condition (0) = 0. From 
the parasupersymmetry relation (3.2b), the Dirichlet boundary condition for Wi implies the Robin 
boundary conditions for Wli (A« 2 V^2)(0) = 0. An important point to note is that, if (X2 = —053, the 
boundary conditions (A~Z W2)(0) = and (A J 1/^2) (0) = are compatible with each other, because 
A^ = — A_Z a for any a. With this choice of the parameter, it follows from the parasupersymme- 
try relation between W2 and 1//3 that the Robin boundary condition (A~ 2 y^)(0) = leads to the 
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Dirichlet boundary condition for y/3, 1^3 (0) = 0- This procedure is easily extended to all \j/ n e H n . 
Basically, the Robin and Dirichlet boundary conditions appear in alternating order. The resultant 
parasupersymmetry invariant boundary conditions are summarized as follows: 



Case N even: < 



Case A" odd: 



'« +1 ^)(0) = 

< 


(n 


= 0,2,4,-- 


■,N- 


2), 


^(o) = o 


(n 


= 1,3,5,-- 


■,N- 


1), 


>+ +1 v«)(o) = o 


(n 


= 0,2,4,-- 


■,N- 


3), 


< v^(o) = o 


(n 


= 1,3,5,-- 


■,N- 


2), 


(A" v„)(0) = 


(n 


= N-1), 







(3.7a) 



(3.7b) 



where {oc„}^~^ are tuned to fulfill the constraints a n = —<x n +\ with n = 2,4,6, ■■■ ,N — 2 for N 
even and n = 2,4,6, ■ ■ ■ ,N — 3 for N odd. Hence, for A" even, we have Af/2 independent parameters 
{cti, (#3, GC5, ■ ■ ■ , CCn-i} in the boundary conditions and parasupercharges, while for A" odd we have 
(N + 1 ) /2 independent parameters { Ot\ , GC3 , OC5 , ■ ■ ■ , CCn-2 , &n- 1 } • 



Type D. Let us move on to the sequence that starts from the Dirichlet boundary condition for the 
bosonic sector, Yo (0) = 0. The discussion is almost the same as the previous one. The result is as 
follows: 



Case N even: 



Case N odd: 



(3.8a) 



(3.8b) 



Yn(0)=0 (n = 0,2,4,- 

(A-v„)(0) = (n = 1,3,5, • 

V„(0)=0 (n = 0,2,4, • 

(A-v„)(0) = (n = 1,3,5, • 

Here the parameters {o^,}^"/ must satisfy the constraints a n = — ce„_i, where n = 2,4,6, •■ ■ ,A" — 
2 for A^ even and n = 2,4,6, - ,N — 1 forA^ odd. We thus have Af/2 independent parameters 
{ai,0£3,0£5,--- ,o;a?_i} for A^ even, while for A^ odd we have (N —\)/2 independent parameters 
{OCi,OC3,(X5,--- ,0,^-2} ■ 



,N-2), 
,N-l), 
,N-l), 
,N-2). 



Parasupersymmetry invariant boundary conditions at x = t are given by just replacing the ar- 
gument x = by x = £ in Eqs. (3.7a)-(3.8b). Having two types of boundary conditions at each 
boundary, we see that there are 2 x 2 = 4 distinct sets of boundary conditions consistent with the 
parasupersymmetry; namely, type (R, R), type (D, D), type (R, D) and type (D, R), where type (R, 
D) for example refers to a combination of type R boundary conditions at x = and type D bound- 
ary conditions at x = t. It should be noted that, due to the parameter constraints a n = —cc n +\ and 
a n = — a„_i (n = 2,4,6, ■ ■ ■ ), there appears only one independent parameter a := (Xi in the type (R, 
D) and (D, R) boundary conditions. 



4 Parasuperspectrum 

In this section we solve the Schrodinger equation H\(f{x) = Ey(x) with the parasupersymmetry in- 
variant boundary conditions and then derive the spectrum. We also discuss the spontaneous breaking 
of parasupersymmetry. 

Throughout this section we will concentrate on the case N even. The case Af odd can be similarly 
analyzed. Both cases are pictorially summarized in Figures 2-4. 
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Figure 2: Schematic parasuperspectrum for the type (D, D) boundary conditions. Parasupersymme- 
try is unbroken. 



4.1 Type (D, D) boundary conditions 

Let us first study the spectrum under the type (D, D) boundary conditions. To this end, let us first 
focus on the positive energy eigenstates, which are Af-fold degenerate thanks to parasupersymmetry. 
The normalized energy eigenfunctions are given by 



Wn,v{x) = \l jsm(k v x)e„ 



(n = 0,2, 4, 2), 



(4.1a) 



Wn,v(x) = (-l)" 2 'y^ + ^( a ) [M(cc n )sm(k v x) + k v cos(k v x)]e n (n= 1,3,5, ••• ,N-\), 

(4.1b) 

where k v = vn/£ (v = 1,2,3, • • • ). We note that relative phases are fixed to satisfy the parasupersym- 
metry relations (3.2a)-(3.2c), or, equivalently, \jr n , v (x) = \jXn=\ ( k l + M 2 (a m ))~ 1/2 (Q + ) n Wo,v(x) 
with the constraints a n = —a n -i (n = 2,4,6,- ■ ■ ,N — 2). The positive energy eigenvalues are given 
by 



V7I\ 2 

-) (v = 1,2,3, 



(4.2) 
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Figure 3: Schematic parasuperspectrum for the type (R, R) boundary conditions. Parasupersymme- 
try is unbroken. 



In addition to these states, there are N/2 distinct negative energy eigenstates 



¥n,o(x) = y & 2M{1^_ I exp[M(« n )x]g» (n = l,3,5,--- ,N-\), (4.3) 
with the energy eigenvalues 

E nfi = -M 2 (a n ) {n = 1,3,5,--, N-l). (4.4) 

Notice that these negative energy states are non-degenerate in general without tuning the parame- 
ters {«] , CC3, as, ■ ■ ■ , CCn-i}- Note also that the negative energy eigenstates (4.3) are the zero-modes 
of the differential operators A Un (n = 1,3,5,- •• ,N — 1) and therefore annihilated by the parasu- 
percharges, Q ± Yn,o(x) = 0. Since the ground state(s) is/are given by the lowest negative energy 
eigenstate(s), we see that in the type (D, D) boundary conditions parasupersymmetry is unbroken. 
Figure 2 schematically shows the spectrum for both N even and odd. 
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4.2 Type (R, R) boundary conditions 

Let us next study the spectrum under the type (R, R) boundary conditions. The positive energy 
eigenstates are N-fold degenerate and given by 



i,vW = (-1)^ k 2 +M^{a + i) [ M ( a "+i) sin (M-^vCos(/: v x)]gn (n = 0,2,4,- ■■ ,JV-2), 

(4.5a) 

^,,vW = (-l) a ^y^sin(*vJc)eH = 1,3,5,- - • ,N-l), 

(4.5b) 

where k v = vn/l (v = 1,2,3, ••• ). Again relative phases are fixed by the parasupersymmetry rela- 
tions with the constraints a n = — a n+ \ (n = 2,4,6, ■ ■ ■ ,N — 2). The positive energy eigenvalues are 
the same as the previous ones 

Ev = (t) 2 ( y = 1 ' 2 ' 3 '---)- ( 4 - 6 > 

Just as in the case of type (D, D) boundary conditions, there are N/2 distinct non-degenerate nega- 
tive energy eigenstates 



Vn,o(x) = y - ^=^^y exp [-M(« n+1 )x] e n (n = 0,2,4, ■ ■ ■ ,N- 2), (4.7) 

whose energy eigenvalues are given by 

E nfl = -M 2 (a n+l ) (» = 0,2,4, 2). (4.8) 

These negative energy states are annihilated by the parasupercharges such that parasupersymmetry 
is again unbroken. The resultant spectrum is schematically depicted in Figure 3. 

4.3 Type (D, R) and (R, D) boundary conditions 

Let us finally solve the model under the type (D, R) and (R, D) boundary conditions. Since these 
two boundary conditions are related by the 7£-transformation, both lead to the isospectral system. 
In what follows we will concentrate on the type (D, R) case. 

Let us first focus on the case M(a)£ > — 1 (or a* := — 2arccot(Mo-C) < a < n). In this case all the 
energy eigenvalues are non-negative and Af-fold degenerate. The normalized energy eigenfunctions 
are given by 



i,vW = (-l) ! 



sin(£ v x)4 (n = 0,2,4,--- ,N-2), (4.9a) 



sin (Jfc v (jc -£))e n (n = 1,3,5, ••• ,N - 1), (4.9b) 



where k v > (v = 0, 1,2, • • • ) is the (v + l)th positive root of the transcendental equation 

M(a) = -kcot(k£). (4.10) 
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Figure 4: Schematic parasuperspectrum for the type (D, R) and (R, D) boundary conditions. Para- 
supersymmetry is spontaneously broken. 



The energy eigenvalues are given by E v = k 2 , (V = 0, 1 , 2, • • ■ ). 

Let us next consider the case M(a)£ < —1 (or — n < a < a*). As shown in Figure 5, in this 
case the ground states energy Eq becomes negative, while all the excited state energies are remained 
to be positive. The ground state wavefunctions are given by replacing ko to iK in (4.9a) and (4.9b) 



Vn,o( x ) = (- 1 ) 

¥n,0(x) = (-1) 



M{a) 



■sinh(?cc)e„ 



(n = 0,2,4,-" ,N-2), 



n-l 
2 



\ K 2 -M 2 (a) l 



smh(ic(x — £))e n (n = 1,3,5, ■ ■ ,iV— 1), 



(4.11a) 



(4.11b) 



where K > satisfies the transcendental equation 

M(a) = - Kcoth(Ki). 



(4.12) 



In this case the ground state energy Eq is given by Eq = — K 2 > — M 2 (a). Opposed to the previous 
subsections, the ground states are not annihilated by the parasupercharges Q . Hence, in the type 
(D, R) and (R, D) boundary conditions parasupersymmetry is spontaneously broken. Figure 4 shows 
a schematic spectrum valid for both N even and odd. 
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Figure 5: Parameter dependence of the energy eigenvalues {Eq,E\,E2, ■ ■■ } for the type (D, R) 
and (R, D) boundary conditions. The blue curve E = E(a) is given by the inverse of the phase 
a(E) = — 2arctan(^cot(\/2?^)) that follows from Eq. (4.10). Spectrum has the period 2k. The 
ground state energy Eq crosses zero at a = a* = — 2arccot(Mo£). 



5 Conclusions and discussions 



It has been long appreciated that supersymmetry and its extensions may be hidden and play a crucial 
role in quantum mechanics with spectral degeneracy. In this paper we showed that the RSKT parasu- 
persymmetries are hidden behind the degenerate spectra in non-relativistic quantum mechanics for 
a single free particle on the carambola graph C#. We imposed the generic cyclic symmetry on 
the graph, which plays the role of grading operator (g^/ivy^ and obtained the Z^-graded Hilbert 
space. We explicitly constructed the parasupercharges and showed that the RSKTparasupersymme- 
try of order p = N — 1 is hidden behind the A^-fold degenerate spectrum. We classified the boundary 
conditions invariant under parasupersymmetry transformations and found that there are only four 
types of parasupersymmetry invariant boundary conditions, of which two lead to the parasuper- 
symmetry breaking. It is interesting to point out here that, in the type (D, R) and (R, D) boundary 
conditions, in which parasupersymmetry is spontaneously broken, the energy spectrum exhibits spi- 
ral structure as shown in Figure 5. Consequently, as one varies the parameter a € S l along S l and 
completes a cycle, the energy eigenvalue E v gets shifted by unit level, E v (a + 2k) = E v+ i(a); in 
other words, nontrivial Cheon's spiral holonomy [16] appears in these cases. Since all the energy 
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levels are Af-fold degenerate, the Wilczek-Zee holonomy [17] might also appear under the adiabatic 
change of parameters that parameterize the basis vectors {e n }. It would be interesting to investigate 
combined phenomena of Cheon's spiral holonomy and the Wilczek-Zee holonomy in this context. 

Before closing this paper it should be mentioned about the relation between our results and the 
previous analysis on hidden supersymmetry structures in quantum mechanics with point interac- 
tions [18-20]. In Ref. [19] Nagasawa et al. considered quantum mechanics on S 1 with two U(2) 
family of point interactions and studied the hidden JY = 2 supersymmetry in the spectrum. Our re- 
sults include theirs when N = 2. In Ref. [18] Uchino and Tsutsui studied quantum mechanics on an 
interval with a single U (2) family of point interactions and revealed the hidden JV = 1 and JV = 2 
supersymmetries in the spectrum. Our results coincide with their JV = 2 results when N = 2 with 
suitable choice of the basis vectors {?o,?i}- They further extended to a system on a pair of two 
intervals each having a single U (2) family of point interactions and showed that large varieties of 
JV = 2 and JV = 4 supersymmetries are hidden in the spectrum [20]. Basically, our N = 4 results 
coincide with their JV = 4 results when the zero-mode energy eigenvalues are tuned to be degener- 
ate. The reason why our results on Cn contain the results on these topologically distinct spaces may 
be understood by topology change [21-23]. It is also interesting to investigate (para)supersymmetry 
structures in topologically distinct graphs in a more systematic way. 
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